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Abstract The baryon content of the universe is known from Big Bang nucleosynthesis and 
cosmic microwave background considerations, yet at low redshift, only about one-tenth of these 
baryons lie in galaxies or the hot gas seen in galaxy clusters and groups. Models posit that 
these "missing baryons" are in gaseous form in overdense filaments that connect the much 
denser virialized groups and clusters. About 30% of the baryons are cool (<10 5 K) and are 
detected in Lya absorption studies, but about half the mass is predicted to lie in the 10 5 -10 7 
K regime, where detection is very challenging. Material has been detected in the 2-5x10 K 
range through OVI absorption studies, indicating that this gas accounts for about 7% of the 
baryons. Hotter gas (0.5-2 xlO 6 K) has been detected at zero redshift by OV1I and OV1II Ka 
absorption at X-ray energies. However, this appears to be correlated with the Galactic soft X-ray 
background, so it is probably Galactic Halo gas, rather than a cosmologically significant Local 
Group medium. There are no compelling detections of the intergalactic hot gas (0.5-10xl0 6 K) 
either in absorption or in emission. Early claims of intergalactic X-ray absorption lines have 
not been confirmed, but this is consistent with theoretical models, which predited equivalent 
widths below current detection thresholds. There have been many investigations for emission 
from this gas, within and beyond the virial radius of clusters, but the positive signals for this 
soft emission are largely artifacts of background subtraction and field-flattening. We discuss the 
various techniques that can be used to detect the missing baryons and show that it should be 
detectable with moderate improvements in sensitivity. 
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1 Introduction and Motivation of the Issues 

It is largely in the past decade that we have come to an appreciation of the missing 
baryon problem and the likely existence of a cosmologically important Warm- 
Hot Intergalactic Medium (WHIM) at low redshift. The belief had been that 
the baryons at high redshift underwent gravitational collapse over cosmological 
time to become the galaxies present today. There was evidence of baryons not 
included in galaxies, such as the X-ray emitting gas in galaxy clusters, where the 
gas usually represented more baryons than the galaxies, but rich galaxy clusters 
were rare and galaxy groups did not possess a bright hot intergalactic medium. 
Also, the Lya forest so prominent at high redshift was sparse near z = 0, so its 
contribution to baryons seemed small. 

A stunning ch ange occurred when a careful census of the baryons at low redshift 
was u ndertaken ( Persic &: Salucci . 19921 ; Bristow Phillippsl . [l99l : lFukugita et al 



Il998h . Thebaryonic mass of galaxies appeared to be only one-tenth of the bary 
onic content at high reds hift, a result first indica ted by studies at least a decade 



earlier (see discussion in iPersic fc Saluccil |1992| ) . The mass at high redshift is 



measured from primordial nucleosynthesis and near z = 3 through the study of 
the Lya lines, which measures the neutral gas content. This baryon content is 
confirmed to high accuracy with recent CMB studies, such a s from the Wilkinson 
Microwave Anisotropy Probe (WMAP; Spergel et al. 20071 ). so it must be low 



redshift baryons that are missing rather than the high redshift baryons that have 
been overestimated. In quoting quantities, we use a flat Friedmann-Lemaitre 
ACDM universe with H = 70 km s" 1 Mpc" 1 , with ft A = 0.72, DM = 0.23, and 
^baryon = 0.045. For this universe, the mean baryon density is 4.2xl0~ 31 gm 
cm -3 , and for an ionized metal- free plasma, the mean particle density is 4.x 10 
cm -3 . The time since the Big Bang is 13.8 Gyr and the distance to a source at 
z = 1 is 3.35 Gpc. 

In trying to account for the missing baryons, there were careful studies of the 
mass in galaxies, and in gaseous form in galaxy clusters, groups, and in the inter- 
galactic medium. Since the first efforts in the late 1990s, there have been many 
observational c ampaigns to improve the numbers. We briefly discuss the up- 
dated census of iFukugita fc Peeblesl (I2004T ). which is consistent with other efforts 



(IShulll . l2003h . ^or the census of the stellar content of galaxies, the availability of 
the SDSS sample was of considerable help as it led to an improved understanding 
of the fraction of s tars i n spheroids compared to disks and to their M/L values 



(iKauffmann et al.l . l200.i l: a large number of low luminosity galaxies is not a vi- 



able possibility for explaining the missing baryons. In estimating the stellar mass, 
there are a variety of models that need to be adopted, such as the initial mass 
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function. iFukueita fc Peebles! (|2004l ) assigns uncertainties to these choices, lead- 
ing to uncertainties of about 25% for the stellar content. Including white dwarfs, 
neutron stars, black holes, and substellar objects, along with the stars burning 
nuclear fuel, this component is about 6% of the total baryon rest mass. In quan- 
tifying the neutral gas content, there have been recent surveys that have greatly 
improved the accuracy of these numbers, suc h as the blind HI P arkes All— Sky 
Survey (HIPASS) that includes 1000 galaxies ( Zwaan et al. . 20031 ). Such surveys 
show that the mass in neutral gas constitutes 1.7% of the baryon content, raising 
the total baryon content in galaxies to 7.7%. More than 90% of the baryons lie 
outside of galaxies, and this became known as the "missing" baryon problem. 

Rich clusters of galaxies have several times more matter in gas than in galaxies. 
The X-ray emission from this hot gas has been used to determine its mass, and 
the baryon content of these rich clusters is consistent with the cosmological value 
(|Allen et all 120021 ) . There are no missing baryons in rich clusters. However, only 



a few percent of the mass of the universe lies in rich clusters, so even though 
gas is abundant, it contributes only about 4% of the total baryon content. For 
decreasing cluster richness, down to galaxy groups, the ratio of the gas mass to 
the galaxy mass declines so that it is less than the mass of the galaxies. The 
difficulty with measuring the gaseous mass in rich clusters or in galaxy groups 
is that the baryon counting stops when the X-ray surface brightness falls below 
detectable levels, at about half of a virial radius for clusters and at smaller radii 
for most groups. In these outer p arts of clu s ters and groups, the implied densit y 
is falling as r _p , where p = 1.5-2 ( Mulchaey . 200d : Rasmussen &: Ponmanl . 20041 ) . 
so the gaseous mass is still increasing at least linearly with radius at the last 
observable radius. Given this uncertainty, it is possible that a significant fraction 
of baryons lie in the outer parts of clusters and groups where the X-ray emission is 
undetectable. While the search for such material is a goal of the field, at present, 
the sum of the galaxies and detected hot gas constitutes only about 12% of the 
baryons. 

One of the ways of measuring the baryon content at higher redshift is through 
the Lya absorbers, a technique that has been applied at lower redshift as well. A 
recent study at z < 0.07 finds that the Lyq-absor bing mat erial accounts for 29% 
± 4% of the baryon content ( Penton et al. 20041 . also see ISembach et"afl 12004 ) . 
While this more than doubles the known baryon content, it still leaves 60% of 
the baryons "missing". These baryons must be difficult to detect, and there 
were a variety of suggestions for the dominant state of the material. However, 
convergence to a prevailing model occurred rather quickly where it was argued 
that gas occupied low overdensity regions in the universe and has a temperature 
m the 10 5 - 10 7 K range. There have been many observations trying to detect 
this hot gas, with some secure measurements, but others being at the detection 
threshold where the reality of the detections have been questioned. 

In this review, we begin by examining the model predictions, where there have 
been significant improvements in recent years. Before reviewing the observations, 
many of which occur at X-ray energies, a brief review of the relevant atomic 
physics is given, followed by a discussion of the X-ray instrumentation, and most 
importantly, the limitations of such observations. This is of special importance 
in regard to the debates of the reality of various X-ray detections. We give a 
critical discussion of each of the major approaches for which the detection of the 
WHIM has been claimed and conclude with the expectations for the future. 
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2 Models for the Warm-Hot-Intergalactic-Medium 



The most accurate models for the large-scale structure formation in a ACDM 
universe have been those that follow only the dark matter, as this is the dominant 
mass component. A gaseous component must be added to model the WHIM, 
along with all relevant heating and cooling processes. Following the epoch of 
reionization near z ~ 6, most of the intergalactic medium has been photoionized to 
a temperature of about 10 4 K. Subsequently, the predominant heating mechanism 
is through the shocks that develop as large-scale density waves collapse in the 
dark matter. Large-scale structure becomes more pronounced with cosmological 
time, so the gas is increasingly shock-heated, reaching temperatures above 10 5 K 
in much of the volume for z < 1. For overdensities of p gas / (/o g as) > 160, virialized 
systems such as clusters of galaxies develop with temperatures of 10 7 -10 8 K. These 
are easily studied by X-ray telescopes, where the curr ent detectability limit is fo r 
Pgas/ (y°gas) > 500, about half of the virial radius (e.g., Pointecouteau et al. 2004! ) . 

The distribution of e clS clS £1 function of temperature and density shows that the 
gas at 10 5 -10 7 K generally lies in regions where the gas overdensity is 10°-10 2,5 , 
although with a fairly wide range, as seen in Figure 1 (|Cen k OstrikeiLl200d ). 
These regions are generally collapsed but not virialized. In addition, there is 
hot gas in very low density regions, which are due to shocks that propagate 
effectively into such environments. Such calculations, even without the effects of 
star formation and stellar feedback, demonstrate t he fraction of gas that would be 
in the 10 5 -10 7 K range (40% in the calculations of Cen &, Ostriker 20061 ; see also 
Croft et al.l l200ll : iDave et al.1 l200ll ) . Unfortunately, most of the observations of 
the WHIM involve non-primordial elements, such as oxygen, so it is necessary to 
include the formation and dispersion of the heavier elements. This necessitates 
having a model for galaxy and star formation, for the supernova rate, for the 
leakage of radiation from the galaxy, as well as for the strength and extent of 
a galactic wind. There is no cosmological simulation that can calculate these 
physical effects from first principles, so they have to be put in by hand with 
educated guesses. 

In their simulations with such feedback, Cen &; Ostriker ( 20061 ) show that the 
energetics of their galactic superwinds is an order of magnitude smaller than shock 
heating by cosmological structures. With galactic superwinds, the predicted mass 
fraction of the WHIM rises from 40% to 50% at z = 0, although the distribution as 
a function of overdensity is nearly unchanged (Figure 2). The extent of a galactic 
superwind establishes a lengthscale over which the metals are distributed, which 
was a value that had to be adopted in prior calculations. This feedback from 
galaxies increases the metallicity of the gas that constitutes the WHIM by about 
two orders of magnitude; their median metalli city is 0.18 Z^, with a broad range. 
A similar simulation, but at higher redshift (jOppenheimer fc Pavel . l200d ) finds 
similar results and they point out that it is difficult to constrain the galaxy wind 
parameters at this point. 

The most common heavy element is oxygen, which has important lines for OVI, 
OVII, and OVIII, so the determination of the various ionic fractions is crucial. 
However, at the very low densities of the WHIM, the recombination times are 
comparable to the Hubble time for nu ~ 10 -6 cm -3 , so non-equilibrium ionization 
effects need to be included, a consideration that s everal researchers h ave examined 
(|Kang et all 120051 : lYoshikawa & Sasaki 120061 : ICen fc EaneL l200d ). There are 
two timescale of relevance: the recombination timescale; and the ion-electron 
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equilibration timescale, which is in usually shorter. Nonequilibrium ionization 
effects will extend the temperature range over which the ions are common. Since 
the WHIM occurs over a wide temperature range, the total predicted absorption 
column of an ion does not change greatly compared to collisional equilibrium 
ionization conditions, but the ratio of the lines from adjacent ionization states 
can change considerably. Therefore, using the OVII/OVIII absorption line ratio 
will not serve as a good temperature indicator, as it does under equilibrium 
conditions. Emission lines are less affected by non-equilibrium effects because 
most readily detectable emission lines occur in relatively denser environments 
(nH > 10~ 5 cm -3 ) where ionization equilibrium is more likely to be achieved. 

There are several goals in making comparisons between the models and the 
data that one hopes can be achieved. The first is the identification of the missing 
baryons, and as a function of temperature, a central prediction (Figure 3). Of 
considerable importance is the location and topology of the WHIM. The models 
predict a density structure, the "cosmic web" (Figure 4), which are the mod- 
erately high ove rdensity regions connecting the easily observed galaxy clusters 
and groups (e.g., Viel et al. 20051 ). Obtaining an observational description of the 



structure of the cosmic web, should it exist, would give us great insight into the 
formation of the WHIM and provide powerful constraints on the models. If struc- 
tural information is available, the densities of the medium can be calculated, also 
of significant interest. In addition, the inclusion of galactic feedback, one of the 
least well-constrained parts of the theory, makes predictions for the metallicity 
distribution of the WHIM, so it would be valuable to confront the models with 
the data. As we will show, observers are struggling with the basic detection of 
the WHIM because the field is very early in its development. 



3 Atomic Physics 

In the UV region, the search for intergalactic baryons has largely focused on 
the detection of the OVI and HI Lya lines. This has been widely discussed, so 
we will just point out a few relevant items. Hydrogen is largely ionized in the 
WHIM, but due to its abundance, it may be detectable at temperatures of up to 
about 3xl0 5 K if it is in collisional equilibrium (ionization fraction of 10 6 ) and 
for total gas columns exceeding 10 19 cm -2 (Figure 5). Unfortunately, this line 
will be broad (FWHM ~ 120 km s _1 from thermal broadening), so detecting it 
upon an AGN continuum, which is often structured, can be a challenge (discussed 
further below). The OVI line reaches a peak fractional abundance near 3xl0 5 
K, and although the fractional abundance does not exceed 0.22, the cross section 
for the 1032 A, 1038 A ground state doublet is large. Near 3xl0 5 K, this line 
is detectable for a total gas density > 10 18 (Z/0.1 Z Q ) cm~ 2 , making it the best 
probe of gas in this temperature range. Beyond 5xl0 5 K, these lines are of 
little use and the X-ray lines become of prime importance. Due to absorption 
by neutral gas in the Milky Way, there is a natural divide between the UV and 
X-ray regions, the former ending at 13.6 eV and the latter becoming important 
above 200 eV. 

For gas at temperatures that emit in the X-ray band, hydrogen and helium 
are fully ionized, so when no metals are present, the only source of emission 
and absorption is by the free- free process. This process is optically thin at all 
X-ray energies, even when calculating the opacity over cosmological distances. 
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The addition of metals greatly increases the emissivity of the plasma since most 
common species have at least one bound electron in the 0.7-10xl0 6 K range. 
Metals contribute to the emission spectra mainly through their bound-bound 
lines which are so plentiful that they form a pseudocontinuum at the typical 
poor resolution of X-ray imaging devices. The strongest of the lines appear 
distinct, with the most important ones coming from oxygen and iron. For the 
study of gas near 10 6 K, the strongest lines are from OVII, a helium-like species 
that has a ground-state triplet at 21.60 A (a permitted line; ls-2p; 574 eV), 21.80 
A (semi- forbidden), and 22.10 A (forbidden). Although only the resonance line 
is important for absorption, the other lines can be more important in emission 
observations. Another important species (at higher temperature) is hydrogenic 
oxygen (OVIII), for which the equivalent of the Ly alpha line occurs at 18.97 A 
(654 eV). Other absorption lines of note are NVII A 24.78, NelX A 13.45, and the 
Ly/? line of OVII at 18.63 A (th is is a highly incomplete list), while Fe XVII and 
related ions are responsible for a variety of lines that are more useful in emission 
observations. 

Even though the product of the peak ionization fraction and the oscillator 
strength is higher for OVII than for OVI, the column needed to achieve a de- 
tectable equivalent width is significantly larger for OVII. This is because the 
fractional equivalent width of a transparent line is inversely proportional to the 
energy of the transition, and as the transition energy for the OVII A 21.60 A line 
is a factor of 50 greater than that of OVI A 1032, the fractional equivalent width 
is significantly smaller. Furthermore, optical instruments have an advantage in 
that they resolve the OVI A 1032 line while the OVII A 21.60 is not resolved by 
current X-ray instruments, typically by a factor of 2-6 (the S/N is proportional to 
the resolving power for an unresolved line). At 0.1 of the solar metallicity and for 
a temperature of 10 6 K, the column needed for a detection toward a strong X-ray 
source is 1-I0xl0 19 cm" 2 . The width of the line is FWHM = 54 (T/10 6 ) 1/2 km 
s _1 if it is at the Doppler width for oxygen, but the sound speed in the gas is 150 
(T/10 6 ) 1 / 2 km s _1 , so if there is sonic turbulence, the FWHM could be closer 
to 200 km s _1 . Also, large scale motion in filaments can easily produce velocity 
widths of a few hundred km s _1 . These line widths imply that the resolution 
needed to resolve the lines are E/(AE) = 1000-3000. However, detections can 
be obtained with a resolution of only a few hundred (the resolution for Chandra 
and XMM near the OVII line is about 400). 

4 X-Ray Instrumentation 

Unlike optical and UV detectors, where one integrates over time and collects many 
photons before readout, X-ray solid-state detectors are used as single-photon 
devices. This is possible due to the relatively slower count rate but there is 
an enormous advantage to the technique in that the number of holes created in 
the CCD is proportional to the X-ray photon energy. Therefore, the energy of 
the input photon is determined to 50-200 eV, leading to a resolution of 10-60, 
depending on the detector and the energy region (0.5-10 keV). This yields a low 
resolution spectrum at every location and permits many narrow-band images 
to be constructed, if desired. When the incoming signal is dispersed through 
reflection or transmission gratings, higher resolution is reached, although with a 
loss of an order of magnitude or more in throughput. At energies below 1 keV, 
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resolutions of 40 0-1200 are achieved with XMM and Chandra (jJansen et al.l . 12001 
Weisskopil2~003h . which is sufficient to begin resolving the lines. 



The sensitivity of these X-ray detectors is not uniform over their energy band 
and the calibration of the detectors can change in unpredictable ways after launch. 
At energies below 1 keV, there is usually a significant low sensitivity region from 
0.28 keV to about 0.4 keV, due to photoelectric absorption from carbon in the 
optical path. At energies below about 0.3 keV, the Chandra Backside Illuminated 
CCDs have no sensitivity, while the XMM Epic camera can detect photons to 
below 0.2 keV, but with declining sensitivity. The changes in sensitivity (the 
calibration) is often monotonic in energy with the lower energies being most 
greatly affected. This is probably due to deposition onto the surfaces of the 
mirror, gratings, or entrance window resulting from outgasing. Recalibration of 
the instrument in flight is a difficult process because there is no ideal celestial 
calibration source, but through clever approaches, moderately good recalibration 
has been accomplished. The importance of proper calibration is that it can lead 
to mistakes in determining the presence of weak emission features, such as from 
the WHIM. 

Another serious issue is the flat-fielding of the data, which must be properly 
accomplished if one is to measure faint extended emission structures. The fields 
of view are typically from 8' to 30', depending on the instrument used. The 
outer parts of these fields are vignetted because the nested mirror design leads to 
one mirror shell shadowing another. This vignetting function can be measured 
accurately for photons coming down the optical axis, but there is a background 
component that is not focused and it can change with time. In addition, there are 
several focused backgrounds: one due to plasma processes in the Solar System; 
one due to diffuse gas in the Local Bubble and above the disk of the Galaxy; 
and one due to AGNs. The first one varies in time and can produce dramatic 
background flar es while the second y aries around the Galaxy, and while that is 
roughly known ( Snowden et all 1997 ). it is not known with precision on the scale 
of Chandra or XMM fields. Separating these components so that one truly has 
a flattened image of the field at energies below 1 keV is challenging and can be 
the limiting step in identifying diffuse structures. 



5 Detection of the 1-5 xlO 5 K WHIM 



At ultraviolet wavelengths, the fractional equivalent width of absorption lines are 
a factor of 30 larger than at X-ray energies (for the same oscillator strengths). 
Thus, the resolution of detectors is often sufficient to resolve the lines (e.g., FUSE, 
and a telescope such as HST has considerably more collecting area than XMM 
or Chandra. These are substantial advantages, and there are two strong lines 
that can be used to trace the WHIM at temperatures below 10 6 K, OVI, and 
the hydrogen Lya line. Under collisional ionization conditions, the lithium-like 
species OVI never becomes the dominant ion, having a peak ionization fraction 
of about 0.25. Furthermore, it is present in only a narrow temperature range (an 
ionization fraction above 0.025 for 2-5xl0 5 K), whereas OVII has a significant 
ionization fraction over nearly an order of magnitude in temperature. However, 
with the ability to detect OVI columns down to 10 13 cm -2 , compared to typical 
X-ray limits of nearly 10 15 cm -2 , and the prediction that there are significant 
amounts of the WHIM at ~10 5,5 K, the UV should be a better hunting ground 
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for detecting the WHIM and studies have borne out this expectation. 

Although OVI had been observed in absorptio n toward hig h er red shift systems, 
the first low redshift detection was ^b^nedbv lTripp et all (|2000h . with a half- 
dozen systems reported by ISavage et all (j2002h . There have been subsequent 
studie s and a large statistical sample to date is presented in IDanforth fc Shulll 
(120051 ). who used FUSE and searched for OVI in 129 known low redshift Lycc 
absorption line systems (z < 0.15) in sightlines toward 31 AGN. They detected 
OVI absorption in 40 systems with column densities of logN(OVI) = 13.0 - 14.5 
and they derive the equivalent width distribution to show that the lower column 
density systems contribute about the same baryonic mass as at higher column 
densities. They determine the space density of OVI (to a limit of 30 mA) of 
dN(OVI) /dz ~ 17 (20% uncertainty), from which they determine the contribution 
to the WHIM of W him = 0.0022±0.003 (Z o /0.1 Z©)" 1 (/ovi/0.2)" 1 , or about 
5% of the baryon content. This is a lower limit to the WHIM from 10 5 ' 5 K gas 
since lower column density absorption lines could contribute signific antly. 



A s lightly larger and somewhat different sample was compiled bv lTripp et al 



(|2006h by using STIS on HST to observe 16 AGNs where they find 44 inter- 
vening OVI absorbers. As the OVI line was redshifted into the STIS bandpass, 
the redshift of the abso r bers (0.12 < z < 0.57) is larger than the sample of 
Danforth & Shulll (120051 ). iTripp et al.1 (j2006h find a slightly larger OVI absorp- 



tion frequency, dN(OVI)/dz = 23±4 for absorption lines with equivalent widths 
greater than 30 mA. Their measured OVI contribution to the baryon fraction is 
correspondingly somewhat greater at 7%. The fraction of WHIM material is pre- 
dicted to decrease with redshift, so the absorption frequency of the set observed 
with STIS might be expected to be 20-30% lower than the set observed with 
FUSE, whereas the opposite is observed. If the model predictions are correct, 
then the absorption frequency dN(OVI)/dz differs at the 2a level between the 
two data sets. It is important to resolve this po ssible discrepancy and a lso to 
extend the surveys to lower equivalent widths, as IDanforth Shulll (120051 ) finds 
that these may add a significant contribution to the mass fraction of gas in this 
temperature range. 

The size of these OVI absorption line samples permit investigators to determine 
whether the material is associated with galaxies. The early studies suggested that 
the OVI lay within a few hundred kpc of a galaxy, but for this statement to be 
made more qua ntitatively, th e chan ce coincidence with random galaxies must be 
accounted for. Stocke et al. ( 20061 ) used the FUSE sample along with several 
galaxy catalogs and they found that for galaxies of luminosity L* or brighter, 
their median separation was 1.8 Mpc, but the distance from a OVI absorber to a 
galaxy was only 0.6 Mpc. Nearly all of the OVI absorbers are found at distances 
of < 0.8 Mpc of a galaxy with L > L* and this gives a maximum characteristic 
sphere of influence around the galaxy for the deposition of the metal-rich WHIM 
near 10 5 ' 5 K. This is a maximum sizescale because smaller galaxies can enrich 
the intergalactic medium, and when galaxies with L < 0.1L* are considered, the 
distance from an OVI absorber to a galaxy falls to 0.34 Mpc. They argue that to 
reproduce the observed dN(OVI)/dz, yet fainter galaxies must contribute to the 
absorption. The sphere of influence of galaxies is consistent with that predicted 
from models. 

Another approach to detecting the WHIM at lower temperatures is through 
neutral hydrogen, and although its ionization fraction is decreasing rapidly near 
10 5 K, because of its much larger abundance, it has a larger column density 
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than OVI for T < 2xl0 5 K (for Zovi = 0-1 Z ). An advantage to using hy- 
drogen is that this is the primary mass component, but the disadvantage is that 
the ionization correction is quite large and the determination of the temperature 
is uncertain. The temperature could be determined from the line width mea- 
surement if thermal broadening were the only broadening mechanism, and there 
are lines observed that are sufficiently bro ad to have temperatures above 10 5 K 
(iRichter et all 120041 : ISembach et all 120041 ): T = lxlO 5 (6/40 km s" 1 ) 2 K, where 
b is th e usual line w i dth p arameter (FWHM = 1.6656). The largest survey is 
that of Lehner et al. ( 20071 ). who estimate that approximately 20% of the local 
baryons may be associated with highly ionized hydrogen. The primary difficulty 
with this method is that the line width can be broadened by bulk or turbulent 
motions in a much cooler gas. Also, multiple unresolved 10 4 K absorption com- 
ponents can blend tog e ther to mimic a broader line from an apparently hotter 
medium. Lehner et al. ( 20071 ) addresses these problems through simulations, and 
while these may resolve such issues, we note that it is notoriously difficult for 
models to properly reproduce turbulence, which can be responsible for the broad 
lines. If the identification of intrisically broad lines can be shown to be truly re- 
liable, the HI Lya absorption would be a valuable tool in measuring the WHIM 
near 10 5 K. 

Finally, an ionization line with a characteristic temperature greater than OVI 
has been detected in the ultraviolet region, NeVIII, although only at the 3.9a level 
(|Savaee et all 120051 ). The NeVIII ion is lithium-like and never even achieves 
a maximum ion fraction in excess of 15%, which occurs near 7xl0 5 K, about 
twice the temperature at which OVI has its maximum ionization fraction. Be- 
cause the ionization fraction distribution for lithium-like species is narrow, there 
is only a very limited temperature range when both OVI and NeVIII may be 
present (at logT = 5.70, /NeVin = 0.021 and /ovi = 0.026 in collisional ioniza- 
tio n equilibrium ) . So while it is possible to produce both lines in the same gas, 
as ISavage etaD (120051 ) suggest, it is equally likely to occur from different tem- 



perature regions where the ioni zation fractio n s are an order of magnitude higher. 
The total column density that ISavage et al.1 (120051 ) suggest would account for a 
NeVIII detection is nearly 10 20 cm -2 , which would be an unusually high-column 
density filament (from theoretical modeling), large enough to produce detectable 
X-ray absorption lines in OVII. This Ne VIII A770.4 line requires a redshift of 
0.3-0.5 to appear in the FUSE band and z > 0.53 to lie in the HST band, so if 
it can be detected at high significance and in other sightlines, it will be a useful 
indicator of the WHIM at intermediate redshifts. 



6 The WHIM Content in the Local Group 



The mass of hot gas in galaxy groups is not well determined because at the last 
radius where X-ray surfac e brightness is u sually measured, the gas mass is rising 
as r m where m« 1.5 (e.g.. lMulchaevll2000h . Therefore, it is possible that galaxy 
groups may possess a cosmologically significant baryon content, and since the 
Milky Way lies in a typical galaxy group, it offers an ideal opportunity for study, 
both through emission and absorption. 

X-Ray absorption line studies are practical for the bright est AGNs and the firs t 
identification of absorption line gas at z « was reported bv lNicastro et al.l( 2002 % ) 
for Chandra observations of PKS 2155-304, shortly followed by lRasmussen et al. 
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(|2003l ). who reported on XMM observations toward 3C 273, Mrk 421, and PKS 
2155-304. Both groups agree that the OVII Ka (21.60 A) and the OV HI Ka 



(18.97 A) lines are detected in PKS 2155-304, althou gh iNicastro et al 



report the detection of a Ne IX Ka (13.45 A) while iRasmussen et al 



(2002) 



report an upper limit. For the OVII and OVIII lines in PKS 2155-304, both 
groups agree that the ion columns are about 5xl0 15 cm -2 in the low opacity 
limit, so the difficulty is how one converts this into a mass of gas. The conversion 
requires values for the ionic fractions and a path length, and it is the choice of the 
path length that can give rise to vastly different values for the inferred gaseous 
mass. 

The first analysis of the X-ray absorption, by Nicastro et al. ( 20021 ). led them 
to conclude that the hot gas fills the Local Group and has a mass of ~10 12 
Mq, significantly more mass than the sum of the galaxies. However, this result 
is model-dependent. They assume that the X-ray absorption lines from OVII, 
OVIII as well as the UV absorption lines from OVI originate in the same gas 
(along with Ne IX). In the X-ray waveband, the lines are unresolved and the 
velocity uncertainty is about 250 km s _1 for the OVII Ka line and 450 km s" 1 
for the OVIII Ka line. At UV wavelengths, the FUSE measurements reveal two 
components at 36 km s _1 and -135 km s" 1 (± 10 km s _1 ), so from the current 
information, it cannot be shown that any of the lines are at the same velocity. 
They consider single-temperatures models of a gas at a variety of densities and 
they include photoionization effects from the X-ray background. They can ac- 
commodate the observed lines in a low density plasma (10 -6 cm -3 ), but the path 
length is 15 Mpc (0.3 Z abundances), which is excluded because a path length 
of this magnitude would show the effects of Hubble expansion, shifting the lines 
in velocity space. To obtain an acceptable solution, they assume that only the 
weaker and broader OVI line (25% of the total column) is associated with the 
X-ray absorbing gas, which permits them to obtain an acceptable solution for n e 
= 6xl0" 6 cm" 3 , Ne/O = 2, and a path length of 3 (0.3/[O/H]) Mpc. This is 
larger than the virial radius for a group of galaxies (1 Mpc) and their path length 
would become 10 Mpc for [O/H] = 0.1, which may be more representative of the 
WHIM. They admit the possibility that there could be a multi-temperature gas 
and that would modify their results s ignificantly. 

A different approach was taken by Rasmussen et al. ( 20031 ) who assume that 
the OVII Ka and OVIII Ka lines are cospatial but that the OVI A 1035 dou- 
blet is from a spatially different region. Unlike INicastro et al.1 (120021 ) . they do 
not detect the NelX Ka in PKS 2155-304 (which had produced a larger Ne IX 
column than either the OVII or OVIII column), but they detect it in Mrk 421, 
with a column (1.5xl0 15 cm -2 ) that is half to one-third of the oxygen columns 
(4.8xl0 15 cm" 2 for OVII and 2.4xl0 15 cm" 2 for OVIII). This reduction in the 
column ratio of Ne/O also reduces the elemental ratio to near-Solar values. The 
observed OVIII/OVII line ratio is used to set the temperature in a collisional 
ionization equilibrium model. They establish a maximum lengthscale of 2 Mpc 
by requiring that Hubble flow not shift the line center beyond its observed value 
(when errors are included). This implies that n e > 5xl0~ 5 cm" 3 and at such 
densities, the ambient X-ray background will not significantly modify the ioniza- 
tion distribution. They argue for a minimum lengthscale of 140 kpc by requiring 
that the soft diffuse emission not be exceeded. They favor a size for the absorbing 
region in the 100-1000 kpc range, which would make this more of a Local Group 
medium than a Galactic Halo medium, with a mass of 10 9 - 10 11 Mq. 
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While both approaches require that their models not exceed the observed dif- 
fuse X-ray background emission, they do not make use of the emission, for which 
a single high resolution spectrum was obtained with a quantum microcalq r imete r 



toward l,b = 90°, 60° and with a 1 sr field of view (iMcCammon et all 1 20021 ). 
This spectrum showed that the strongest single feature was the OVII Ka triplet 
and this yields an emission measure for OVII, given a temperature and metal- 
licity. When combining the absorption and emission measures for OVII and for 
an ion fraction of 50%, the length scale becomes about 20 kpc, with n P = 9xl0~ 4 
cm" 3 , or a gas mass of 4x10 s M (jSanders et all 12002 : Bre gman <fe Llovd-Davied . 
2007). A similar result is obtai ned by using th e emis sion measure in ferred from 



the R osat broad-band data of Snowden et al. (2000), as shown by Fang et al 
(2006). In this analysis, the gas is in a Galactic halo and the gas mass is an order 
of magnitude less than the neutral gas content of the Galaxy. This gas is most 
likely to have come from the disk, heated by supernovae, so the assumption of 
Solar abundance for the metallicities is justified. For this model to be correct, 
some, if not most of the OVIII absorption comes from a slightly hotter gas (0.1- 
0.2 dex), but most realistic systems have temperature fluctuations. Also, if the 
filling factor is less than unity, the halo would be larger. A minimum size to the 
halo can be derived by assuming an ionization fraction of unity for the OVII, 
which leads to a lengthscale of 5 kpc. 

Subsequent studies have sought to cl arify which models are most likely. The 
largest lengthscale model advo cated by Nicastro et alj ( 20021 ) is unlikely to be 
valid, based on several studies. IWilliams et al. (2005) shows that, for a sightline 
toward Mrk 421, the OVI A1032 absorption line has a larger Doppler parameter 
than for the OVII A21.60 line. They further argue that this Doppler parameter 
is inconsistent with the OVI(HVC )/OVII ratio (they m ake similar arguments for 



an observation to ward Mrk 27 9; Williams et al. 2006). For the same sight line 



toward Mrk 421, Savage et al.1 ( 20051 ) shows that the OVI absorption between 
-140 km s _1 and 60 km s" 1 is associated with low-ionization gas in the Galactic 
thick disk/halo. For the higher velocity gas, from 60-165 km s _1 , absorption by 
CIII is also seen, and based on the column density ratio N(OVI)/N(CIII), they 
argue that this gas cannot coexist with OVII and OVIII at the same temperature. 

Some researchers argue that the OVII and OVIII absorbing gas is produced 
in a Galactic Halo rather than a Local Group medium. They analyze a set 
of z absorption lines detected toward a sample of the brightest AGNs us- 



20051 : iFang et all 120061 : 



ing both Chandra and XMM o bservations (I Wang et al.1 _ 

Bregman fc Llovd-DaviesI 20071 ). One of the distinguishing properties is whether 
the absorbing gas has a redshift similar to the Milky Way (0 km s _1 ) or to the Lo- 



cal Group (-250 km s 1 ). This difference is less than the instrumental resolution 
(HWHM « 400 km s _1 ), but a line center can be determined more accurately 
than the instrume ntal resolution. For a Chandra observation toward LMC X-3, 
IWang et al.1 (|2005T ) finds a velocity of 56 km s 1 (-83, +194 km s 1 ), more s ugges - 
tive of a Galactic origin, but with large errors. Bregman &; Llovd-DaviesI ( 20071 ) 
used the data from the four highest S/N AGNs observed with XMM , where the 
random uncertainty in the calibration of a single observation is 111 km s _1 . They 
find a mean velocity of 70 ± 55 km s -1 , which is consistent with the Milky Way 
velocity but 4.6<r different from the Local Group velocity. Th e velocities found 
by ot hers are consistent with a velocity near km s _1 (e.g., iRasmussen et al. 
20061 1. 

Another argument is given by ( Fang et al. . 20061 ) who show that if the absorp- 
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tion lengthscale were 1 Mpc, and if this were typical of galaxy groups, intergalac- 
tic absorption would easily and commonly be detected as the sightlines of AGNs 
pass through relatively common galaxy groups. Their constraint indicates that 
the lengthscale be less than about 0.5 Mpc (at the 99% confidence limit) and this 
might be improved by making this test with the higher S/N XMM data. 

Additional insight may be gained by determining whether the distribution of 
absorption line strengths has a Galactic dependence or a Local Group depen- 
dence. If there were a halo around the Galaxy for which the halo radius is <J 30 
kpc, lines of sight across the Galaxy would have greater absorption than those 
away from the Galaxy. The Local Group has a different signature caused by 
its elongation toward the MW-M31 axis, with the greatest column toward M31 
and the least in the anti-M31 direction. By good fortune, the anti-M31 direc- 
tion is toward the bulge of the Milky Way (l,b = 301°, 22°), where the Galactic 
Halo model predicts the greatest absorption columns. This makes the two models 
nearly ortho gonal, which is the ideal si t uation . To define the angular distribution 
of columns, Bregman &: Llovd-Daviesl ( 20071 ) used a sample of 26 extragalactic 
lines of sight, 17 of which have uncertainties in the OVII Ka equivalent width 
of 10 mA or less (the median equivalent width is about 20 mA). The data show 
no enhancement in the sight lines within 30° of M31, nor is there a minimum in 
the equivalent widths toward the anti-M31 direction. However, there is a cor- 
relation between the equivalent widths and the Galactic 3/4 keV background, 
which measures the emission from gas in and around the disk, including a halo 
(Figure 6). This adds to the evidence that the OVII gas is primarily of Galactic 
origin, probably of radius 10-50 kpc. A hot halo of similar size is seen in the 
recent Chandra X-ray observation of the massive spiral NGC 5746, where X-ray 
emission extends to at least 20 kpc. The estimated gas mass is 3xl0 9 M Q for 
a metallicity of 0.04 Zq. This metallicitiy is uncertain and as the derived mass 
depends approximately inversely with the metallicity, the mass would be several 
times smaller for a metallicity of 0.2 Z . This mass estimate is similar to, but 
larger than the value inferred for the extended halo around the Milky Way. 

A final consideration that bears upon this is observations involving the Magel- 
lanic Stream, cold gas that has been tor n from the LMC/SMC system by the grav- 
itatio nal tidal forces of the Milky Way ( Mastropietro et al. . 20051 : Connors et al 
2006). There are indications that the Magellanic Stream is interacting with a 
dilute medium, both from ionization and structural considerations. The Stream 
is composed of a variety of clouds and the surfaces of the clouds in the pre- 
dicted direction of motion (along th e Stre am toward the LMC/SMC) are bright 
in Ha. This led Weiner &: Williams ( 19961 ) to argue that the emission was shock 
heated due to an interactio n with a gas where n > 10 -4 cm -3 (more recently 
and throughly discussed by Putman et al.1 120031 ) . If t he Stream clouds are pres - 
sure confined by the hot medium, n <^ 3xl0~ 4 cm" 3 (Stanimirovic et aD . 20021 ) . 



Lower estimates to the density are given by ISembach et al. (2003) who argues 
that some of the OVI absorption detected by FUSE occurs at the boundaries 



cm 



between clou ds (T <; 10 4 K) and the 10 6 K medium (with n ^ 10 4 -10 
Also, iMuralil (|2000l ) suggests a density < 1 cm if the clouds are to survive 
evaporation by thermal conduction, but as Mastropietro et al. ( 20051 ) point out, 
the density limit may be an order of magnitude higher for a lower (but sensible) 
choice of temp erature. Also, in their gr avitational and hydrodynamic modeling 
of the Stream, iMastronietro et, al.l «) argue that a model with hot dilute gas 
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improves the results, where they choose a mean density of 8.5 xlO -5 cm -3 within 
50 kpc. A total hot gas density of 10 -4 cm -3 would lead to an electron column 
density of lxl0~ 19 cm~ 2 , which would only contribute 20% to the mean inferred 
OVII Ka equivalent width for Solar abundance (and a mass o f 6x10 s Mm). 



At greater radii, there is a useful density estimate from iBlitz fc Robishaw 



(2000) who point out that the nearer dwarf galaxies are devoid of cold gas while 



the more distant ones sometimes have HI. They attribute this difference to ram 
pressure stripping and derive a value for the gas density of 2.5 xlO -5 cm -3 at 
a distance of 200 kpc. This also corresponds to an electron column density of 
lxlO -19 cm -2 , so it would never dominate the OVII Ka equivalent width, al- 
though the gas mass would be lxlO 10 M . This gas mass is about 15% of the 
baryonic content within 200 kpc, so it is an important but not dominant compo- 
nent. 

To summarize, these studies indicate that the Milky Way has a hot gaseous 
halo that produces the observed OVII Ka absorption and the gas has a with a 
temperature near 10 6 K and a mass of 4xl0 8 (//20 kpc) 2 M . A more dilute 
gaseous medium probably extends through the Local Group and with a mass of 
~10 10 Mq, which provides only a modest contribution to the total baryon content. 
There are fewer constraints on the gas containing OVIII, as its absorption line 
is less commonly detected. It may arise in the a Local Group medium or in 
the hotter regions surrounding the Galaxy, if there are regions where the gas 
temperature is near 10 63 K. 

7 X-Ray Absorption Detection by the WHIM? 

The detection of the WHIM in the X-ray lines would reveal gas at 10 6 -10 7 K and 
open a new line of research that is complementary to the detection of cooler gas 
through the UV lines. This is a pioneering area that has been the focus of much 
attention and observing time, which is justified by the magnitude of the issue 
involved. Our belief is that the evidence for a discovery of this importance must 
be compelling and by this criteria, it is unlikely that there has been a detection 
of dilute 10 6 -10 7 K beyond the Local Group. 

The primary claim for the detection of the 10 6 -10 7 K WHIM comes from the 
observation of Mrk 421, the brightest AGN at X-ray wavelengths, and since this 
is a featureless BL Lac object, it is an ideal target against which to search for 
absorption lines. As discussed above, the strongest line is expected to be the 21.6 
A resonance line of OVII, although other lines, such as the 18.97 A Lya line of 
O VIII could be strong for temperatures near 10 6 -10 7 K. 

iNicastro et al. ( 20051 ) reported on two long (96-100 ksec) Chandra observations 



of Mrk 421, one taken with the LETG in combination with the High Resolu- 
tion Camera (HRC) and one with the LETG in combination with the ACIS 
CCD array. These observations were obtained when Mrk 421 was particularly 
bright, so they obt ained 2.7xl0 6 photons (HRC) and 4.3xl0 6 photons (ACIS; 
Kaastra et al]l2006h . or a combined amount of about 5500 photons per 50 mA bin 



near 21.6 A. They combined the data from the two instruments for their analysis. 
Ideally, one would like to fit a simple spectral model, such as a power-law with 
Galactic abso rption, and then iden tify spectral features in the normalized spec- 
trum. When Nicastro et al. ( 20051 ) tried this, enormous residuals result. They 



show that these residuals are due to calibration errors near elemental absorption 
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edges. To deal with this calibration problem, they fit a power law plus absorption, 
but the abundances of the absorbing gas were permitted to have values signifi- 
cantly different from those of Galactic interstellar gas. This procedure is meant 
to account for contaminating material in the optical path of the instrument in 
addition to the usual absorption by Galactic interstellar gas. A much better fit 
was obtained with C and Ne abundances about three times the Solar value while 
N and O were nearly absent. Following this procedure, there were still a number 
of significant deviations in their spectral fit, so they added 10 broad Gaussian 
features, six in emission and four in absorption. This fitting procedure, which 
includes unphysical abundances and arbitrary Gaussian components is used to 
produce a flattened spectrum. From this flattened spectrum, they begin a search 
for physically meaningful absorption lines. 

This flattened spectrum requires that additional line components be added, 
because the quality of the fit was unacceptable, having a \ 2 of 2892 for 1598 
degrees of freedom. They point out that the largest remaining deviations were 
clearly identified with absorption at z = 0. Consequently, they added 24 absorp- 
tion line fits and this reduced the \ 2 to 1911 for 1529 degrees of freedom; this 
is the final normalized data set they work from. Although this is a significant 
reduction in x 2 > it still is not an acceptable fit, with a probability of occurrence 
much smaller than 10~ 4 . 

They examine whether any of the 24 absorption lines lie at redshifts between 
that of Mrk 421 (z = 0.030) and z = and they find two weak absorption line 
systems at redshifts of 0.011 and 0.027. The strongest line is the OVII 21.60 A line 
at z = 0.011, for which their conservative significance is 3.8a. This significance 
results from a Gaussian fit where all Gaussian parameters and the continuum 
are free to vary. There are no other lines at this redshift above the 3a level, 
but they identify two nitrogen lines (NVI Ka and NVII Ka) at the 3.0cx and 
3.1a level and the OVII 21.60 A line (2.8a) with a redshift of 0.027. An equally 
strong absorption line at 24.97 A has no identification. The OVI A1032 line is not 
detected at these redshifts with FUSE. This is an important non-confirmation 
because the OVI line should exist at about the same temperature as the NVI 
line and FUSE is sensitive to columns two orders of magnitude smaller than this 
Chandra observation. A HI Lya line is detected at z = 0.01016, but the column 
is orders of magnitude smaller than that needed to produce the X-ray N or O 
lines, although this is not necessarily a conflict as the ionization temperatures of 
these species ar e quite different . Bas ed on the column density and frequency of 
the X-ray lines, iNicastro et al. (|2005h calculate a cosmological mass density for 
the WHIM that would be consistent with model predictions. If this is correct, it 

constitutes the discovery of the hot missing baryons. 

The detection of these lines has been questioned ( Kaastra et al. . 20061 ; Rasmussen et al. 



2006), mainly using XMM data but also using 244 ksec of Chandra calibration 



data for Mkn 421 taken about a year after those by Nicastro et al. ( 20051 ). The 
sum of the calibration data provide a spectrum of 3.2 million photons, about 
20% more than the LETG/HRC spectrum and 2 5% fewer photons than the 
LETG/ACIS spectrum reported bv INicastro et al.1 (|2005l ). iKaastra et all (|2006h 
reanalyze the Chandra data, including these calibration observations and an ad- 
ditional LETG/HRC observation from 2000. Instead of fitting a spectral model, 
they fit a 12-point spline with nodes separated by 0.5 A, plus they removed the 
strongest instrumental features and Galactic absorption. They tre ated separately 



the LETG/HRC data set, the LETG/ACIS data obtained bv INicastro et al. 



Bregman 



15 



( 20051 ). and the LETG/A CIS calibration obser vations. At the location of the 
strongest line reported by iNicastro et al" ( 20051 ). the OVII 21.60 A line at z = 
0.011, they confirm a dip in the spectrum in two of the three data sets. This fea- 
ture was absent in the LETG/ACIS calibration observations (Figure 7). They fit 
the instrumental line spread function to the residuals of the normalized spectrum, 
rather than fitting a Gaussian. T hey find that only the LETG/ACIS spectrum 
obtained by lNicastro et al. (|2005h shows a detection above the 2a level (at 2.5c"), 
an d only for the OVH li ne at z = 0.011. Their fit to both sets of data discussed 
by Nicastro et al. ( 20051 ) is 2.7<r, lower than the 3.8<r quoted by them, with the 
presumed difference being due to the somewhat different fitting procedures. 

Processing the XMM data poses several additional challenges because there are 
a number of features in the detector, such as hot pixels, which are not addressed 
when using the standard SAS processing. By using any recent version of the SAS 
(2006 or later), the resulting spectrum has several instrumental features seen in 
every bright spectrum. One of these happens to lie near 21.83 A, which is very 
close to the OVII 21.60 A line at z = 0.011 (21.85 A), thus mak i ng it impossible 
to de termine its absorption line strength ( Ravasio et all 20051 ; Williams et al. . 
2006). However, if one constructs a complete map of all problematic pixels and 



uses that to define the acceptable data, it is possible to obtain high-quality spec- 
tra over nearly the enti re XMM-RGS range, and that is the procedure used by 
Rasmussen et al. (120061 ) . That effort describes the data processing and analysis 



of 955 ksec of RGS data from mid-2000 through late-2005, and by good for- 
tune, Mrk 421 was quite bright during many of these observations. This led to 
a net spectrum where there were about 26,000 counts per 50 m A bin, about 
five times more than th at used by lNicastro et al. I (l2005h in their Chandra data. 
Rasmussen et al. ( 20061 ) fit their RGS data similarly to the Chandra data and for 
the OVII line at z = 0.011. They find no detection. When they form a weighted 
average that includes the Chandra and XMM data, the line is a l.Ocr feature, 
w ell below the thresho ld for a detection. 

Kaastra et al.1 <|2006h argue that the statistical significance of the lNicastro et al 



(2005) results should be because they would have accepted an absorption at 
any wavelength over a redshift range of 0.03. As many instrumental resolution 
elements fit in this range, the statistical significe must include an appropriate 
number of trials, which is not considered in th e discovery paper. We find the 
reanalysis and criticism of kaastra et al J « to be convincing and conclude 
that there is no statistically significant detection of intergalactic absorption in 
the OVII line toward Mrk 421. The rms of the XMM data is smaller than for 
the Chandra data, and the derived upper limit is logN(OVII) = 14.6. 

There is another claimed detection of absorption by the WHIM toward H1821+643. 
This is about an order of magnitude fainter in soft X-rays than Mrk 421, but it is 
at z = 0.297, so there is more red shift range to search. It was observed for 470 ksec 
with the Chandra LETG/ACIS ( Mathur et al. . 20031 ) and the resulting spectrum 
had 2xl0 5 photons, about 50 times less than the number of photons accumulated 
for Mrk 421 from all LETG observations. The observers searched for X-ray ab- 
sorption that lay within 0.025 A of the redshifts of six intervening OVI absorption 
line systems determined from UV studies. They searched for absorption by OVII 
Ka, OVIII Ka, and NelX Ka, and they report one system associated with OVI 
absorption with a S/N > 2 (at 2.2), and two systems with S/N > 2 that lie 
within ±1000 km s _1 of the OVI system (one in OVII Ka and the other in OVIII 
Ka). Given the ranges used by these investigators, there were effectively nine 
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line resolution elements to find matches with the OVI systems and 20 trials for 
the matches near OVI systems. For the number of effective trials, finding a few 
2(7 absorption features has at least a 10% probability of oc currence accor d ing to 
our calculations. The more detailed error simulations of Kaastra et al. tod ) 



would suggest that the significance is rather poor. These claimed features to- 
ward H1821+643 are not convincing detections eit her and along with th e failure 
of XMM to confirm the earlier Chandra result of Nicastro et al. I (l2005h in Mrk 
421 leads us to conclude that the 10 6 -10 7 K WHIM has not yet been detected 
beyond the Local Group. 

In this ongoing search for the hot WHIM, there has been some confusion re- 
garding the capabilities of the Chandra LETG compared to the XMM RGS. For 
example, it is often mentioned that the LETG has significantly higher resolution 
than the RGS, which quickly compensates for the difference in collecting area. 
This is a true statement at some wavelengths, because for a spectrally unresolved 
line where the continuum is well-defined, the S/N of a detection is linearly pro- 
portional to the spectral resolution. However, it increases more slowly with the 
effective area, being proportional to its square root. At 21.6 A, the spectral reso- 
lutions are nearly the same for both observatories, being 0.05 A (FWHM) for the 
Chandra LETG and 0.06 A for the XMM RGS. Also, the wings of the line spread 
function contain more flux for the XMM RGS than for the Chandra LETG. In 
contrast, XMM has an effective area that is about a factor of 3-4 greater, which 
more than compensates for the modest spectral resolution advantage of Chan- 
dra. For the particular case of Mrk 421, t he XMM observations are about five 
times longer tnan the Chandra data used hv lNicastro et alJ «). Although the 
brightness differences during their observations approximately compensated for 
the collecting area difference of the two telescopes, the longer exposure time with 
XMM led to a spectrum with a lower rms. 

The failure t o detect the OVH l ine in the WHIM is not a problem for cosmo- 
logical models (|Cen fc Fand . 120061 ). The frequency of OVII absorption columns 
exceeding 10 15 cm -2 is about dN/dz ~ 1.5 (Figure 8), so a Az range of about 1 
is needed for a reasonable chance of making a detection and Mrk 421 only offered 
A z = 0.03. 



8 X-Ray Emission From The WHIM in and Near Galaxy Clus- 
ters 

There has been much activity in the past decade in the effort to identify diffuse 
emission from the missing baryons both in and outside of galaxy clusters. The 
emission projected within the virial radius of galaxy clusters is discussed sepa- 
rately from the emission beyond the virial radius. Any emitting material must 
have a fairly high overdensity (>10 2 ), so even if it is truly detected, it is not the 
typical WHIM material predicted by models, but the upper density tail. 



8.1 Soft Excess Emission Within the Virial Radius of Galaxy 
Clusters 

Clusters of galaxies are among the brightest X-ray sources, containing hot gas 



10 7 -10 8 K) that accounts for more baryons than the visible galaxies (e.g. jAllen et al 



2002). This gas has been studied at X-ray wavelengths for a variety of reasons, 



such as to determine the baryonic and gravitating mass or to study abundance 
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patterns. The instruments used are spectroscopic imaging devices (proportional 
counters or CCDs) where each incoming photon can be assigned an energy to 
modest accuracy (3-50%, depending on the device and the photon energy). This 
permits one to carry out spatially-resolved spectroscopy, whereby a model spec- 
trum is fit to the spectral energy distribution in some region, such as an annulus 
around the cluster center; the data typical cover the energy range 0.2-10 keV. 
The spectrum is dominated by free-free radiation for typical cluster energies (3- 
10 keV), with line emission being of secondary importance. Another important 
parameter of the fit is the absorption by cool Galactic gas (< 10 5 K), as well as 
by material within the emitting source. This absorbing gas is due to the sum of 
the HI, H2, and warm ionized gas at ~10 4 K, and as only the HI can be measured 
accurately, the absorption column is usually a fitted quantity, although with con- 
straints. For typical column densities of 10 20 -10 21 ' 5 cm -2 , this absorption occurs 
at the low-energy part of the X-ray spectral energy distribution, 0.1-0.7 keV. 
Therefore, the simplest spectrum has a single temperature, abundance values, 
and an absorption column density. If the simplest spectrum produces an unac- 
ceptable fit, more complicated fits are adopted. Unfortunately there are several 
technical challenges, especially at the lowest energies, where instrumental cali- 
bration is most difficult and where background subtraction can pose problems. 

The results of spectral fitting for galaxy clusters had led to claims of both excess 
absorption and excess emission. The claims about exc ess absorption beg an with 
Einstein Observatory Solid State Spectrograph spectra ( White et al.l . ll99ll ). where 



the soft emission was less than would be expected from Galactic absorption and 
a single-temperature cluster spectrum. This implied a significant amount of gas 
cooler than 10 5 K within the cluster. However, these SSS spectra had to be 
corrected for the buildup of ice in the optical path. There was concern that if 
the correction for ice was wrong, the modified spectrum would have a soft X-ray 
deficit. This so ft X-ray deficit was no t confi rmed with subsequ ent instruments 
such as Rosat (jArabadjis Bregmanl . 2000) or XMM-Newton (jPeterson et al" " 



2003), so we can safely conclude that the original study was incorrect and that 



there is no substantial absorbing medium within clusters. 

Wheterh there is an additional emission component at soft X-ray energies (0.1- 
1 keV) has been a much more persistent and controversial subject. This line of 
study began with the detection of the Virgo Cluster by the Extreme Ultraviolet 
Explorer (EUVE), an instrument that was not designed to detect extragalactic 
objects. Howev er, one of the imaging de vices on the EUVE is the Deep Survey 
(DS) telescope ( Bowver fc Malina . 1 99 ll ) . It is sensitive to photons up to about 
0.19 keV, the low energy part of the soft X— ray region. At these energies, the 
optical depth is significant (r=3 at 0.16 keV), but still small enough so that a 
modest amount of X-ray emission can pass through the gaseous Galactic disk 
and be detected by the EUVE. For a Galactic column density of 2xl0 20 cm -2 , 
typical of high Galactic latitude sightlines, the effective bandpass of the EUVE 
i s 0.144—0.186 keV (FWHM) but with a collecting area of only about 0.5 cm 2 
( Bregman et al. . 20031 ). iLieu et ah (1996) took the existing spectral fits to the 
X-ray data and asked whether the model also explained the EUVE measurement, 
but found that the model underpredicted the EUVE observations. They claimed 
that much of the emission detected by the EUVE was due to a new and previously 
unknown component and this was the beginning of soft excess emission from 
galaxy clusters. 

If this soft excess emission is due to nonthermal emission from cosmic rays in 
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the galaxy clusters ( Sarazin Lieu . 19981 ). it has no impact on the presence of 
a WHIM. However, if this emission is thermal, it must be due to gas at 1-3 xlO 6 
K, and the mass of gas involve d must be compara ble to the hotter ambient 
gaseous mass in the cluster (e.g., Kaastra et al. 20031 ). Therefore, it would be of 
cosmological importance. 

The claim of a soft excess detected with the EUVE was reexa mined by the 
mission P.I., S. Bowyer, and collaborators. iBerghofer et al. (|2000h argued that 
flat-fielding correct i ons w ere not properly applied to the EUVE data (see also 
Bowver et al. 19991 . l200ll ). They constructed flat fields from several very long 



observations and compared them to the same clusters for which excess emission 
had been claimed. They found that the flux detected by the DS on EUVE could 
be explained by a model that fit the cluster X-ray emission, with no additional 
soft c omponent; t he on ly exception was the Coma cluster. In response to this 
work, Lieu et al. ( 19991 ) argued that there could be problems with using such a 
blank field, since time-d ependent changes in the detector background would lead 
to analysis errors. Also, iLieu et ah I (11999T ) suggested that the choice of threshold 
cutoffs in the pule-height data co uld lead to inco r rect r esults. These criticisms 
where considered and dispelled by Berghofer et al. (2000), who found no evidence 
for time variation in the instrument's background and no significant effect due 
to the choices made in the pulse-height threshold level. To summarize, with the 
exception of the Coma cluster, Bowyer and collaborators do not confirm the soft 
excess that Lieu and collaborators commonly detected. 

Although it is remarkable that the EUVE was able to detect some galaxy 
clusters, it has far less collecting area than either Rosat or XMM in the same 
energy regions. The effective collecting area of the DN on EUVE is about an 
order of magnitude less than the collecting area of the Rosat PSPC or of the XMM 
EPIC pn (t hin filter or open posi tion) for the same energy band and absorption 
conditions (IBreg man et all liB . Therefore, most studies of the soft emission 
have concentrated on the soft X-ray observations of clusters. Some studies by Lieu 
and collaborators, using the Rosat data showed that they were consistent with the 
EUVE data in that some fraction of the soft emission (E < 0.5 keV) required a 



new soft component in addition to the hotter cluster component (|Lieu et all Il996l : 



Kaastra et al. . 19991 : Bonamente et al. . 200ll ). These researchers also showed that 



the radial surface brightness profile of the soft excess was flatter than the surface 
brightness distribution of the hotter gas. 

However, an examination of the Rosat data by lArabadjis Bregman (Il999h 
found that for annuli outside the central region, the cluster spectra could be fit 
with a single-temperature hot thermal spectrum plus Galactic absorption and 
that no additional soft component was needed. Again, the exception was the 
Coma cluster. The primary difference between this and the work where a soft 
excess was claimed has to do with the treatment of Galac tic absorption. Rather 
than fixing the Galactic absorption at the 21 cm value, lArabadiis &: Bregman 
(1999) permitted the value to vary within the uncertainties of the 21 cm mea- 
surement. The measurement of the 21 cm HI column is not straightforward 
because stray radiation can enter the radio beam and this has the greatest affect 
on sightlines out of the plane where the columns are smallest. It is possible to cor- 
rect for this sidelobe contamination, but the uncertainties in th e corrections and 



in calibrations lead to the uncertainties in the 21 cm column (Hartmann et al 



199d : iMurphv et ail ll996T). Also, it wa s important to use the improved He ab 



sorption cross sections ( Yan et al. . 19981), because He is the primary absorber of 
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soft X-rays at these columns. 

There is likely to be a systematic effect in the absorption of soft X-rays, because 
of fluctuations in the Galactic interstellar medium. This systematic effect can 



cause an apparent soft excess (jBregman et all 120031 ). The Galactic HI column 



is known to have small-scale structure and variations, yet when correcting for 
Galactic X-ray absorption in extended sources, one applyies an average HI column 
measured from 21 cm emission. At low optical depth, it is sufficient to equate 
the X-ray absorption column to the mean column measured from 21 cm emission. 
However, at energies where there are moderate and high optical depths, such as 
the 0.15-0.5 keV region where investigators search for the soft excess, the X- 
ray absorption calculated from the mean HI column (as measured from 21 cm 
emission) leads to a systematically incorrect spectral fit. This can be seen if we 
consider two lines of sight with colums N+5N and N-<5N. The average would be N, 
which is what the 21 cm observer records, but the average fractional transmitted 
X-ray flux is l/2(exp(-o"E (N+5N))+exp(-o"E (N-5N))), where <te is the absorption 
cross section at some energy. At low optical depth, the transmitted flux equals 
F Q (1 - cje N), as though there were no column density fluctuations. However, 
when the optical depth exceeds unity, the above expansion is not accurate and 
the full exponentials must be retained. The extra radiation transmitted in the 
lower density line of sight is greater than the extra amount absorbed in the other 
line of sight so that the net transmitted flux is greater than if there were a 
single uniform column (i.e., l/2*(exp(-<TE (N+<5N)) + exp(-cTE (N-5N))) > exp(- 
<7e N)). This effect leads to a correct fit at low optical depths but a positive 
residual at higher optical depths that would be interpreted as excess emission. 
To calculate the magnitude of the effect, one would need an accurate map of the 
actual fluctuations in HI within the aperture being considered, and such data 
are generally not available. Instead, we estimated the effect based on N(HI) 
variations seen in other parts of the sky and find that it can be a contributor to 
the soft excess at moderate optical depths, typically for N(HI) > 3xl0 20 (E/0.3 
keV) 3 cm~ 2 . 

Following several studies with Rosat, the satellite Beppo-SAX was used by 
IKaastra et al. I ljl99fll ) to obtain data for Abell 2199 and to argue that it contains 
a soft component. In a rather forceful criticism of this work, iBerghofer Bowyer 
(120021 ) examined the same data set but used a different approach to their analysis. 
They found no evidence for an additional soft component either in Abell 2199, 
or in Abell 1795, in a g reeme nt with the pri or analysis for the s ame clusters by 
Arabadiis Bregmanl (|2000l ). In a rebuttal, Kaastra et al. ( 20021 ) argue that the 



analysis by Berghofer and Bowyer was oversimplified, leading to their failure to 
detect the soft excess, which they reaffirm. The differences between these works 
once again have to do with the technical details of background subtraction and 
flat-fielding. 

More recently, Chandra and XMM-Newton give us another opportunity to 
examine the problem. Due to its larger field of view and lower energy resp onse , 
XMM proves to be better-suited for this problem, so Kaastra et al. (2003) used 
XMM data to search for the soft X-ray excess emission in 14 galaxy clusters. 
When adopting the Galactic HI column as the X-ray absorbing gas, a few clusters 
show a soft X-ray deficit, which can be understood if there is additional gas along 
the line of sight, which could be due to molecular hydrogen. In the majority of 
the clusters, they find evidence for excess soft emission and they show that it 
is broadly extended across the clusters. Several of the best cases are discussed 
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in more detail and they attribute some of the soft excess to emission from the 
OVII line. In addition, they point out that this line appears to be redshifted 
and is therefore not associated with the Galaxy or Local Group. The presence of 
OVII emission would point to the emission mechanism being thermal and from a 
relatively cool component (0.1 keV) in or around the cluster, such as the WHIM. 



The best cases showing a diffuse excess were examined separately by lBregman fc Lloyd-Davies 



(2006), since they noticed that the strength of the soft excess emission seemed to 



correlate with Galactic latitude. A better correlation was found between the soft 
excess and the soft X-ray background near 1/4 keV (tak en 5° from the cluster ; 
the X-ray background used was the Rosat R12 bands of ISnowden et al. (|l994l ): 



this is in turn correlated with Galactic latitude; Figure 9), so it was suggested 
that the background subtraction method may have led to the soft excesses. The 
reason for this is that most of these clusters are too large for one to use a back- 
ground from the same image. Therefore, they used a standard background that 
had been compiled from a variety of observations around the sky, which represent 
an average X-ray background. Above 1 keV, the X-ray background is dominated 
by AGNs, and this component is fairly isotropic on the sky. However, below 1 
keV, the X-ray background is dominated by the emission from Galactic gas and 
this emission can vary significantly around the sky. If the standard background 
below 0.5 keV is less than the Galactic value toward an individual cluster, then 
removal of the standard background inevitably results in an apparent soft X-ray 
excess. 

To test this explanation, iBregman fc Llovd-Davlesl l|200d ) constructed back- 



grounds appropriate for each galaxy cluster by using backgrounds from other 
fields that had similar values for the Galactic soft background, N(HI), and the 
instrumental particle background (for clusters Abell 1795, Abell 1835, MKW 3s, 
and Abell S1101, also known as Sersic 159-03). When these individual back- 
grounds were subtracted from the observations of clusters, the resulting spectra 
could be fit without the need for an additional soft component. Also, one cluster, 
Abell 1835, is small enough that an on-chip background can be used and the 
resulting spectrum did not show soft excess emission. 

In summary, the reality of soft excess X-ray emission in galaxies has been 
seriously questioned, whether the data come from the EUVE, Rosat, Beppo-SAX , 
or XMM-Newton. In every case, the disagreement can be traced to issues of flat- 
fielding (exposure maps) and background subtraction, which are tricky problems. 
Given the scientific importance of discovering the WHIM in galaxy clusters, we 
believe that the proponents of the soft X-ray excess must give unambiguous 
evidence for its presence. Until this is accomplished, the detection of the WHIM 
in galaxy clusters cannot be claimed. 

8.2 The Coma Cluster 

The Coma cluster is a possible exception in that the cluster spectrum appears 
to require extra emission at low energies (< 0.4 keV) even if one uses a Galactic 
HI column near the low end of its likely range (in this general direction, N(HI) 
= 0.9±0.1xl0 20 cm -2 ). All investigators agree t hat an excess is present in this 
cluster, with the most recent work in X-rays by Finoguenov et al. ( 20031 ). who 



used XMM-Newton data. They use a background from another part of the sky, 
which they scale to the HI column appropriate for the Coma cluster. This has 
the shortcoming that the soft X-ray background varies around the sky, as does 
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the fraction from the Galactic halo and the Local bubble, so there is no unique 
method for scaling the background from one region to another. It is not possible 
to use a local background because the virial radius of Coma (assumed to be 3 
Mpc) is 1.7°, much larger than the field of view of the EPIC camera. Given 
these shortcomings, these researchers have performed a careful analysis and in 
their resulting spectrum, they detect the OVII and OVIII features, from which 
they obtain a redshift of 0.007±0.004±0.0015 (best-fit, statistical, and systematic 
error). This is within 1.3a of zero redshift, expected for Galactic emission, but 
it is about 3a from the redshift of the Coma cluster, 0.023. In support of soft 
emission from Coma, there are weak Ne IX and OV III absorption line s at 2.3cr 
and 1.9a, obtained with the RGS on XMM-Newton (jTakei et all . 120071 ) . If these 



could be confirmed, it would be strong support of cooler material around the 
Coma Cluster. 

There are three possible explanations for this apparent excess. The first is 
that the Galactic soft X-ray background happens to be high in this direction 
because of the properties of Galactic gas. The Galactic HI column lies in a local 
minimum in this direction (very close to the North Galactic Pole), which certainly 
was not caused by the Coma cluster. A minimum in N(HI) will let more Galactic 
Halo emission through, leading to a brightening that is unrelated to the Coma 
cluster. Such brightenings in the R12 (1/4 keV) band are seen in other HI minima 
locations, including the Lockman Hole. This local hole near the Coma cluster can 
be see n in the dust extincti on map (Figure 10) that utilized the IRAS and COBE 
data ( Schlegel et al. . 19981 ). The typical extinction about 4° from the center of 



Coma is about 0.015 magnitudes. Even this is several times smaller than the 
typical sightline out of the Galaxy. However, the Coma cluster lies in a region 
where the extinction is more typically 0.008-0.009, an extremely low value. 

A related possibility is that N(HI) is actually lower than measured because 
the sidelobe contam ination has been underestimated, an issue t hat is worst for 
the lowest columns (jHartmann et al" . 19961 : Murphy et al. . 19961 ). With a mean 
column of 6-7xl0 19 cm -2 , a spectral model does not require a soft excess compo- 
nent, but this value lies about 3a below current measurement for the HI column 
in this direction. 

The third possibility is that there is a soft component and that it is produced 
by inverse Compton scattering of the Cosmic Microwave Background by energetic 
particles in the cluster. Such a possibility is attractive because this is one of the 
only clusters with a bright nonthermal radio halo. This m o del h ad been sug- 
gested previously ( Sarazin &: Lieul . 19981 ). but Bowver et al. ( 2004 ) argues that 
the EUV emission, compared to the X-ray emission, falls more steeply than would 
be expected if the nonthermal electron density decreased a s the hot gas density . 
A steeper decrease in the nonthermal particles is implied, so Bowver et al. ( 20041 ) 
argues that the nonthermal particles responsible for scattering photons are sec- 
ondaries caused by interactions of the cosmic ray primaries with the hot gas. 

To summarize, there appears to be an excess of soft X-ray emission from the 
direction of the Coma cluster, but it may not be from a ~ 10 6 K medium within 
the cluster. The soft X-ray excess could have a Galactic origin or it may be due to 
non-thermal processes within Coma. Higher quality X-ray spectral observations 
could help to clarify the matter. 
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8.3 The Absence of OVI emission 

The soft diffuse emission has been interpreted as being due to thermally emitting 
gas with a temperature in the range 0.5-1 xlO 6 K, which is close to the peak of 
the cooling function. Unless balanced by an unknown heating mechanism at this 
temperature, the gas will continue to cool through the 3xl0 5 K range, where 
most of the thermal energy is emitted through the OVI AA 1032 1038 lines. This 
emission is accessible with certain extreme UV spectros copic instruments, notabl y 
the Hopkins Ultraviolet Telescope (HUT) and FUSE foixon et all Il99fil . l200lh . 



the latter having better sensitivity. Five clusters were observed with HUT: the 
Hercules cluster; Abell 1795; Abell 1367; the Coma cluster; and the Virgo cluster. 
Only upper limits were obtained to the OVI emission, which were inconsistent 
with the presence of 5xl0 5 K gas in the Vir go cluster. Their m ore recent FUSE 
observations, of the Virgo and Coma clusters ( Dixon et al.l . [200ll ) are significantly 



deeper. They also fail to detect the OVI emission. These observations rule out all 
published models in which the soft excess is caused by thermal emission in these 
galaxy clusters. The absence of the OVI emission further brings into question 
the presence of WHIM material in the Coma and Virgo clusters. 

8.4 X-ray Emission Beyond the Virial Radius of Galaxy Clusters 

In principle, a successful emission study of the WHIM could measure the tem- 
perature, intensity, metallicity, and most importantly, the topology of the gas 
(expected to be a web of connected filaments). In practice, this is challenging 
because one is observing a volume emission measure that is expected to be much 
fainter than galaxy clusters or even galaxy groups. For a typical rich cluster where 
X-ray emission is easily detected, a representative density is n e ~ 3xl0 -4 cm -3 , 
or p g as/(Pgas) ~ 1500. The density is falling roughly as r~ 2 , so at double the 
radius the gas density has dropped by a factor of four but the emission measure 
has decreased by more than an order of magnitude and is becoming very difficult 
to detect (for p gas / (p ga , s ) < 500). The typical overdensity in the WHIM is 10°-10 3 
and for gas above 10 K (typical of groups and clusters), it is 10 2 -10 4 . So if one 
were to detect most of the mass of the WHIM, it would be necessary to reach 
overdensities of 10, about 2000 times lower surface brightness than currently de- 
tectable limits. In addition, for clusters or groups, one can obtain a brightness 
in an annulus, improving photon statistics, in contrast to filaments where the 
location, shape and extent are unknown. An additional reduction in the surface 
brightness would occur if the metallicity in the filaments is lower than in galaxy 
clusters, because most of the emission is from metal lines at these temperatures. 

Aside from the issue of the "soft excess" from clusters, where the presence or 
absence is often due to the method used to flat-field the data, there are other 
efforts to detect this emission. 

Emission from the WHIM with temperatures below about 0.1 keV will be 
impossible to detect due to Galactic absorption, so the part of the WHIM that 
one seeks to detect is the 10 6 -10 7 K range. As there is a weak correspondence 
between temperature and overdensity, gas in this temperature regime would have 
a characteristic overdensity of 10 2 (range of 10 1 - 10 3 ). 

There have been several discussions about the emission around and between 
galaxy clusters. Usi ng Rosat data, whi ch provides a larger field of view than any 
current instrument, Soltan et al. ( 19961 ) argued that Abell clusters have extended 
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low surface brightness halos that extend to a characteristic radius of 7 Mpc, about 
twice the virial radius (using the Rosat All-Sky Survey). However. iBriel Henry 
(1995) searched for filaments connecting 54 Abell clusters, with separations of 4- 
17 Mpc. Hierarchical clustering theory predicts that the emission from the WHIM 
should be brightest along loci between clusters, but lBriel &: Henryl ( 19951 ) failed 
to detect emission. For a temperature of 1 keV (1.2xl0 7 K), a metallicity of 0.3 
of the Solar value, and a filament width of 0.5 Mpc, their 3cr upper limit to the 
electron density is 9xl0 -5 (1/0.5 Mpc) -1 cm~ 3 , where I is the assumed thickness 
of the filament. For a temperature of 0.5 keV (6xl0 6 K) and a metallicity of 0.1 
Solar, the density limit would be twice as large. These upper limits still represent 
an overdensity of 400-1000, so the failure to detect such filaments is not in conflict 
with model predictions. 



Sea rches for filaments were made using pointed Rosat PSPC data and lScharf et al 



(|2000h have claimed the detection of a filament at about the 3<r level. Its surface 
brightness is about three times lower than the limit of Briel Henry ( 19951 ). cor- 
responding to an electron density that is about a factor of two lower. There are 
some concerns with this result, acknowledged by the investigators. For example, 
the location of the cluster is superimposed on one of the bri ghtest regions in the 
soft 1/4 keV X-ray background (the R12 band). Although lScharf et al.1 (|2000h 
work in the 0.5-2 keV band, the plasma responsible for the Galactic 1/4 keV 
background also has strong OVII emission at 0.58 keV, so much of this would fall 
in their band. Without more detailed spectral information, it will be difficult to 
determine whether this filament is related to structure in the Galactic soft X-ray 
background or whether it is intrinsic to the CL 1603 field. 

A search for filamentar y connections between c lusters in the Shapley super- 
cluster was undertaken by Kull &: Bohringer ( 19991 ). where Rosat data were used 
to analyze a 6° by 3° region containing three luminous central clusters (Abell 
3562, Abell 3558, and Abell 3556). About 1.5 Mpc between Abell 3558 and 
Abell 3556, dilute emission was detected and it is softer than the bright cluster 
emission, consistent with g as in the 0.5-1 keV ra nge. The surface brightness is 
about twice the 3cr limit of Briel &: Henrvl (jl995l ). This appears to be a strong 
detection. However, the projected position places this within the virial radius of 
both clusters, making it difficult to claim that this is a unvirialized filament of 
moderate overdensity. 

More recently, XMM data were used to investigate an extended 4 Mpc filament 



near Abell 85 that was originally discovered with the Rosat PSPC (jDurret et al 



120031 ) . They confirm the Rosat result and show that the extended feature has the 
same position angle as the major axis of the central cD, the cluster galaxies, and 
the nearby galaxy clusters groups. The temperature of the filament is measured 
to be 2 keV (2.5 xlO 7 K), and although this is cooler than the clusters, it is 
significantly hotter and brighter than expected for the WHIM. Also, its projected 
length extends only a bit beyond the virial radius of Abell 85 (about 3 Mpc), 
so this may be interacting withe the cluster, as the researchers suggest (among 
other possibilities). 

One of the greatest challenges of this sort of work is removing the instrumental 
background, which can be larger than the sought-after signal. In that regard, 
the newest X-ray telescope, Suzaku, is superior in that it has a lower soft X- 
ray instrumental background than previous missions. The spectral resolution of 
the CCD is also improved, allowing narrow wavelength band studies to focus 
on emission from particular ions, such as O VII. Among the first results is a 
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searc h for OVII and OVIII emission around the cluster Abell 2218 (iTakei et all 



20071 1. They examine regions 10-20 Mpc from the cluster but do not detect excess 



emission, with a limiting val ue that is abou t a fac tor of five below previous claims 
of soft excess observations ( Kaastra et all 20031 ). Their density limit is n# < 



7.8 xlO -5 cm -3 for Z = 0.1 Zq and a path length of 2 Mpc. This corresponds to 
an overdensity of 270, which is higher than the typical overdensity of the WHIM, 
so this upper limit is consistent with expectations. 

In conclusion, although there does appear to be faint emission at about the 
virial radius of rich clusters, there is no definitive detection of a filament that is 
well beyond the virial radius and with a temperature in the 10 6 -10 7 K range. 

8.5 Shadowing the Diffuse Component of the X-Ray Background 

Rather than detect faint filaments directly, another approach is to detect the sum 
of many filaments. This has the advantage that the surface brightness will be 
many times greater than a single filament, especially since some filaments will 
be seen along the long axis. There have been calculations as to the expected 



surface brightness and the angular fluctuations around the sky (jVoit et all 120011 ; 



Croft et alll200ll ). The emission is dominated by multiple lines of Fe (L series) 



and O (OVII and OVIII), for abundances above 0.01 of the Solar value. The 
important lines lie below 1 keV and there are two selection effects that lead to 
the emission being confined to the 0.2-1 keV range. Emission below 0.2 keV is 
absorbed by the Milky Way and the Galactic X-ray background is bright below 
about 0.4 keV. Therefore, emission from high redshift gas will be redshifted into 
these unfavorable energy ranges. Also, at higher redshifts, the fraction of WHIM 
gas decreases, leading to less emission. Calculations of the cumulative emission 
shows that it should be most prominent around 0.7 keV with a moderate range, 
so a u seful band i s 0.4- 1.0 keV. The emission-weighted mean redshift is <z> 
0.32 ICen et all (Il995ll and the peak effective temperature is about 3xl0 6 K 



dCroft et all 1200 1). 



In this energy band, three components contribute to the diffuse X-ray back- 
ground: the Galactic soft background; unresolved AGNs; and the WHIM. The 
WHIM contribution is expected to be the minority component , contributing at 



the 1-10% level, although this predic tion is model-dependent (ICen et all Il995l : 
Croft et all l200ll : IPhillips et ail l200lh . A very important non-astronomical com- 



ponent is the background caused by the instrument and by material in the Solar 
System (cosmic ray, scattered Solar X-rays, etc.). Usually, this is comparable 
to the astronomical diffuse background and it varies in time. Separating these 
many components and confidently extracting the small component of the WHIM 
is challenging and can only be accomplished under special circumstances or with 
a specially designed (future) mission. 

One of the special situations in which instrumental backgrounds can be re- 
moved is the Rosat All- Sky Survey, where there are many overlapping and con- 
tinuous fiel ds. By using the an t icorre lation betweenm the RASS and the Galactic 
HI survey, iKuntz &: Snowdenl (j200(j ) were able to separate local emission from 



emission beyond the HI layer. They were able to model the emission beyond 
the HI layer with two thermal components at l.lxlO 6 K and 2.9xl0 6 K. This 
emission contains all sources beyond the disk, including unresolved point sources. 
Based upon deep point source studies, they remove these components and arrive 
at a uniform all-sky value of 7.5 keV s -1 cm~ 2 sr~ _1 keV" 1 in the 3/4 keV band 
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(|Kuntz et all l200ll >. Some portion of this emission is probably associated with 



the Galactic halo, which could be the dominant component. This constitutes the 
upper limit to the 3/4 keV emission from the WHIM. It is comparable to the 
value suggested by ICen Qstriker of 7 ke Y s 1 cm 2 sr 1 keV 1 , but 

an order of magnitude larger than the estimate of Phillips et al.1 ( 200ll ) for the 
WHIM. 

A method of directly determining the extragalactic background, separated from 
other backgrounds is through X-ray shadowing. The technique uses a cloud that 
lies at some interesting distance and for which the opacity due to photoelectric 
absorption is significant enough to be measured. A spatially large diffuse emis- 
sion component behind the cloud will appear less bright at the cloud location 
because of absorption by the cloud. That is, the cloud produces a shadow. This 
was used successfully to demonstrate that a significant amount of the 1/4 keV 
diffuse background lies beyond the Local Bubble. It can also be used to shadow 
the diffuse emission from the WHIM, provided that sufficiently dense clouds lie 
beyond the Local Group. Spiral galaxies have the requisite amount of cold gas 
but they also contain emitting sources (X-ray binaries, supernova remnants, etc.). 
Therefore, it is necessary to find a situation where the emitting sources are not 
unimportant. Such a situation can occur for edge-on galaxies, either because gas 
that extends beyond the disk of stars or because the column density of gas is so 
large in the disk that it absorbs essentially all of the emission from sources in 
the galaxy. This approach is made easier with Chandra as one can identify and 
exclude point sources to good accuracy. 

The ed ge-on galaxy NGC 891 w ould appear to be a suitable object for this 
approach. Bregman &: Irwin ( 20021 ) used the high optical (and X-ray) extinction 
region of the midplane of NGC 891 for shadowing; this region subtends less than 
1 square arcminute. Aside from disk point sources, the galaxy has some halo 
emission and if this is on a scale larger than the HI disk, the shadowing approach 
would not work. They reported a possible detection of a shadow. In the ex- 
tinction region, 68.4 photons were expected but only 51 photons were detected, 
which is a deficit that would be expected to occur only 1.7% of the time. How- 
ever, in a subsequent and significantly longer observation, the shadow was not 
confirmed (Bregman, in preparation; Figure 11). The upper lim it provided by 



this observation is less restrictive than that of Kuntz et al 



(2001). 



9 Other Techniques for Detecting the WHIM 

Although this review has focused on the X-ray and UV absorption lines as the 
primary tools for detecting the WHIM, there are other promising techniques. 
Most will require instrumental improvements to be competitive. However, as 
there are plans for significant instrumental improvements in the future, they 
could become important new approaches for WHIM studies. 



9.1 Dispersion Measures 

An ionized plasma has a frequency-dependent index of refraction, which pro- 
duces a measurable effect at radio wavelengths. The effect is well-known from 
the studies of pulsars, where the pulse at lower frequencies lags the same pulse 
observed at higher frequencies. The difference of the arrival times of the same 
pulse at different frequencies is proportional to the dispersion measure, which is 
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the integral of the electron density over the path length. The electron column 
density is precisely the quantity that one would like to measure through cosmo- 
logical space. From it, one can determine the mass content of the WHIM without 
the abundance or ionization corrections that must be confronted with most other 
methods. The dispersion measure has been measured through the Milky Way 
as well as for pulsars projected upon the Magellanic Cl ouds. The total electron 
column out of the Milky Way is a bout 5xl0 19 cm -2 ( Taylor Corded . Il993l : 



Corded . 2001 : Crawford et all 12001). This is similar to th e column inferred from 



the OVII Ka absorption (iBregman Llovd-Daviel 120071 ) 



The 20 pulsars in the Small and Large Magellanic clouds ha ye dispersion mea- 



sures that are several times the value due to the Milky Way ([Manchester et al 
2006). Also, these values vary by up to a factor of five between pulsars, so both 
the higher values of the dispersion measure and the variation is attributed to the 
ionized interstellar gas in the Magellanic Clouds. The expected dispersion mea- 
sure from diffuse gas between the Milky Way and the Magellanic Clouds is less 
than variation between the Magellanic Cloud pulsars, preventing us from using 
them to gain information about halo or Local Group gas. 

Unfortunately, pulsars have not been detected beyond the LMC, so we cannot 
use this approach to determine the total electron column in the Local Group, 
or beyond. Studies of AGNs have failed to find pulses or other variations on 
a sufficiently short timescale that can be used to determine dispersion measures 
over cosmological distances. Another class of candidates are gamma- ray bursters , 
provided that they have sufficiently sharp prompt radio emission, liokal (120031 ) 



has examined this possibility and discussed that it may be feasible with a future 
instrument such as the Square Kilometer Array. One problem with using either 
AGNs or gamma-ray bursters for dispersion measure determinations is that the 
environment of these objects may produce dispersion measures in excess of the 
cosmological values, severely compromising the technique. 

9.2 Radio Hyperfine Lines 

Another approach may be provided by using the hyperfine lines of highly ionized 
species. These are the equivalent of the 21 cm hyperfine transition for hydrogen 
that occur in either hydrogenic or lithium - like s pecies with significant nuclear 
magnetic moments ( Syunyaev Churazovi . Il984l ). Unfortunately, most of the 



common elements are composed of alpha particles, which have weak magnetic 
moments and do not produce useful hyperfine lines. Nitrogen is one of the only 
exceptions to this rule because it is made in abundance through the CNO process. 
A hyperfine line occurs at 53.2 GHz in hydrogenic nitrogen (NVII), but the 
atmosphere is opaque at this frequency. Consequently, this line cannot be used 
to probe the Local Group, but it can be used at modest redshift (z > 0.1) where 
the atmosphere is transparent. In this line, the emission from the WHIM is 
too smal l to be detectable, but th e detection of this line in absorption holds 
promise ( Goddard &: Ferland . 20031 ). For HI, the hyperfine absorption feature 



is greatly reduced by stimulated emission, since the two levels are close to their 
Boltzmann ratio, with the collision time being much shorter than the lifetime of 
the excited level. However, for NVII, the transition time is less than the time 
between collisions for WHIM conditions, so stimulated emission is unimportant 
and it is a much better absorber. 

The NVII hyperfine line can be compared to the OVII Ka absorption, where 
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a column of 10 16 cm -2 yields a good absorption feature. Nitrogen is six times 
less common, so the central optical depth of the line, which is easily resolved 
with radio receivers, is < 10~ 3 , and for a good detection, it would be necessary 
to achieve opacities of 10 -4 or less. One difficulty is that few continuum sources 
are strong in the 40-50 GHz range, but there several sources brighter than a 
few Jy. We observed several radio-loud AGNs to search fo r such weak lines 
in 3C 273, 3C 279, 3C 345, and 4C 39.25 (jBregman fc Irwinl . l2007l ). While no 



absorption features were discovered above the 5a level, these observations showed 
that it was possible to obtain an rms per channel as low as 10~ 4 with about 
an hour of integration. For integration times of 25-30 hours for the brightest 
objects, the sensitivity will be comparable to long X-ray observations (300 ksec) 
of the brightest sources, provided that the rms of the radio spectrum continues 
to decrease as the inverse square root of the integration time. Radio observations 
have their own set of systematics that must be dealt with, such as standing waves, 
but these systematics are different from the X-ray observations, and thereby 
provide a valuable complement in the search for the WHIM. 



A v ariety of other hyperfine lines are available for study, but as lGoddard <fe Ferland 



(2003) point out, nearly all of the others are too rare to have any likelihood of 



being detectable in absorption by the WHIM. 

9.3 The Sunyaev-Zeldovich Effect and the WHIM 

The passage of the cosmic microwave background (CMB) through a hot cloud 
of electrons will cause the photons to gain energy, leading to a distortion of the 
spectrum, the Sunyaev-Zeldovich effect (S-Z). At frequencies below about 200 
GHz, this leads to an apparent decrease in the temperature of the CMB, and 
the magnitude of the shift is proportional to the line integral of the pressure 
of the hot medium. Th is effect is se e n tow ard rich clusters of galaxies through 
dedicated observations ( Grego et al. . 200ll ). and even in the survey data from 



the first year of Wilkinson Microwav e Anisotrpy Probe (WMAP), it is detected 
( Hernandez-Monteagudo et all 120041 ). 



In practice, there will be a S-Z effect caused by the WHIM in the form of 
filaments, but this will be much fainter than for clusters. Relative to clusters, 
the temperature in the WHIM is at least an order of magnitude lower and the 
column density through a filament is also at least an order of magnitude smaller. 
Therefore, the S-Z effect in WHIM filaments will be at least two orders of magni- 
tude smaller than for rich clusters. Ha nsen et al. I (|2005h searched for such a signal 



by using galaxies to trace WHIM filaments, and correlating the galaxy density 
with WMAP data. They conclude that the S-Z signal from the WHIM cannot 
be detected with the current WMAP data and they derive upper limits. 

Although direct detection of WHIM filaments may be out of reach at the mo- 
ment, it may be possible t o probe the outer parts of clusters and groups, although 
this is challenging as well. Hernandez-Monteagudo et al. ( 2004! ) sear ched for but 



failed to see a S-Z signal from structures less bound than clusters. iMvers et al 



(2004) report on a similar study that used the same WMAP data in which they 



have a marginal detection of the S-Z effect on a size scale (5 Mpc) somewhat 
greater than the virial radius of the cluste r. Given th e disag reement between 



the groups and the marginal detection by IMvers et al.l (|2004l ) , the evidence is 
not compelling that the S-Z effect has been detected beyond the virial radius of 
galaxy clusters. However, future instruments may provide enough improvements 
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in sensitivity to make such studies feasible. 



10 Final Comments and Future Prospects 

The search for the WHIM is a relatively new field in which the majority of studies 
are pioneering efforts that are technically challenging. We can look at the final 
census and try to anticipate where future instrumentation might make the great- 
est impact. The combination of the stars and gas in galaxies, galaxy clusters, and 
galaxy groups accounts for about 12% of the baryons, while he cool gas that pro- 
duces Lya absorption at low redshift accounts for about 30-50% of the baryons. 
Of the remaining material, there is a secure detection of WHIM material in the 
temperature range 1-5 xlO 5 K through the OVI absorption lines seen in a number 
of lines of sight, and this comprises about 7% of the baryons. The remaining 30- 
50% of the baryons are predicted to lie in the temperature range 0.5-30xl0 6 K. 
Although there have been claims for the detection of cosmologically significant 
amounts of gas in this temperature range, these results are not compelling. We 
define a compelling detection as one that is above 5cr, that is confirmed by multi- 
ple instruments (when appropriate), and does not disappear for different choices 
of data reduction procedures (such as using different backgrounds or flattening 
procedures). 

Gas near 7xl0 5 K can be studied if the redshifted resonance lines from Li-like 
Ne VIII A773 can be detected. The equivalent width of these lines are about 
an order of magnitude less than the OVI A1035 lines, so with sufficiently long 
observations with FUSE or HST (with the Cosmic Origins Spectrograph), it 
might be possible to detect this spe cies above the 5<r level (a 3.9a detection is 
currently claimed; Savage et al. 20051 ). Above this temperature, progress requires 



X-ray observations, typically from the resonance lines of hydrogenic or He-like 
species, such as OVII and OVIII. 

At temperatures near 10 6 K, the OVII Ka line is the most promising diagnostic, 
while the OVIII Ka line is most common in gas at 10 6 - 35 K, but with a wide 
temperature range (up to about 5xl0 6 K). These species have been detected 
in absorption by gas at zero redshift. This absorption is most likely caused by 
Galactic halo gas, which is not massive enough to be of cosmological importance. 
There is no compelling evidence for absorption by these ions from intergalactic 
gas. This is not in conflict with models, which show that about an order of 
magnitude greater sensitivity is needed, and for several lines of sight, in order 
to detect such absorption. This improvement in sensitivity is provided by the 
present design of Constellation- X , which should usher in a rich period of discovery 
for WHIM studies. Studies of absorbing gas at 10 K and above require using 
species such as Mg XII (the Ka line), which are at least an order of magnitude 
less common than oxygen and where the absorption line occurs at higher energies. 
Detecting these absorption lines will be challenging. 

We conclude that there is no compelling evidence for emission from WHIM 
material within clusters or beyond the virial radius. The most important ob- 
servations would be detecting the WHIM beyond the virial radius. This would 
reveal the shape and density of the cosmic web of baryons. Any observation 
along these lines must demonstrate that the gas is at the expected redshift, so 
this will require sensitive spectroscopic observations with very low instrumental 
noise and a moderately broad field of view. This cannot be achieved with current 
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instrumentation, but is perfectly feasible with future wide-field instruments using 
quantum microcalorimeters. 

Finally, there are other technologies that may prove useful, such as S-Z mea- 
surements and radio NVII absorption line observations. The S-Z measurements 
are most sensitive to hotter gas, as the signal is caused by the pressure integral 
along ones line of sight, so it would nicely complement X-ray absorption line ob- 
servations. However, at least one order of magnitude improvement in sensitivity 
will be required for this approach to be feasible. 

Uncovering the missing baryons is a feasible goal, requiring only sensitivity 
improvements that can be attained with existing technologies. We hope that 
these sensitivity improvements will be realized in the coming decade through 
the construction of the next generation of instruments. The result would be a 
watershed of new discoveries. 

JNB would like to thank Jon Miller, Jimmy Irwin, Renato Dupke, and Mary 
Putman for many helpful conversations. Also, Andy Rasmussen and Fabrizio 
Nicastro helped to clarify some technical issues, while Ken Sembach and Pat 
Henry provided valuable comments. John Kormendy and Roselyn Lowe- Webb 
had many suggestions that helped to claify the paper in the editing process. JNB 
would like to acknowledge support from NASA through the LTSA grant program, 
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Figure 1: The mass distribtuion of the warm gas (T < 10 5 K), the WHIM (10 5 -10 7 
K ) and the hot gas (T > 10 7 K) function of gas overdensity in the simulation 
of ICen fc Ostrikerl (|20Q6l ). Most of the mass of warm gas lies at overdensities of 
0.03-30, whereas most of the WHIM has overdensities in the range 1-1000, and 
most of the hot gas has overdensities in the range 30-10 4 . 
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Figure 2: The evolution of the different baryonic components of the universe 
from z = 3 to z = 0, from the model of lCen & Ostrikerl tOOd ). where galactic 
superwinds (GSW) are included (the evolution of the WHIM component is also 
shown without GS Ws) . This is typical of other simulations that show the present- 
day WHIM component to be about 50% of all baryons. 
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Figure 3: The differential mass fr action as a function of temperature for the sim- 
ulations of ICen & Ostrikerl (|2006h . where the solid curve includes galactic super- 
winds whereas the dashed line does not. The mass of the WHIM is fairly evenly 
distributed between 10 5 K and 10 7 K in the galactic superwind model. There is 
a local minimum in the mass fraction near 3xl0 5 K, the peak temperature for 
the OVI absorption line diagnostic. 
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Figure 4: Th e density distribution f or 10 5 -10 K gas from a cosmological ACDM 
simulation of lCen fc Ostrikerl (|200fil ). where = 0.37, ft 6 = 0.049, a 8 = 0.80, and 
where L = lOOtr 1 Mpc. Most of the mass of the WHIM lies within the filaments 
that connect the higher density region. 
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Figure 5: The ion fraction distributions, represented as column densities for 
a total gas column of 10 19 cm -2 and metallicities of 0.1 Zq. The UV lines 
are effective at detecting absorbing gas for T < 5xl0 5 K and currently have 
significantly better sensitivity than the X-ray OVII Ka and OVIII Ka lines, which 
are good diagnostics for gas temperatures in the range 0.5-5xl0 6 K. Absorption 
by OVII, OVIII and NelX have been detected at z = 0, probably because of 
the higher metallicity of Galactic halo gas. The NVII line (dotted), which has a 
hyperfine line in the radio region, is also shown. 
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Figure 6: The OVII equivalent width as a function of the Galactic R45 inten- 
sity (3/4 keV) where the 26 sam ple objects have been binned into seven groups 
( Bregman Lloyd-Daviea . I2007I ). The correlation of the two quantities suggest 
that most of the OVII absorption has a Galactic origin, probably caused by a 
halo of size 15-100 kpc. 
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Figure 7: The 21.0-21.6 A s pectrum of M r k 421 obtained with the Chandra 



Observatory and reduced by feaastra et al. ( 20061 ). The top spectrum is the 



LETG/ACIS observation when Mrk 421 was brightest, the middle spectrum is 
the sum of the LETG/HRC-S observations and the bottom spectrum is the more 
recent calibration observations taken with the LETG/ACIS (it is brighter than 
the middle spectrum but has been displaced downward by a constant factor for 
clarity). This spectrum includes the OVII Ka resonance line at 21.60 A, which is 
easily detected at z = 0. The dotted l ine sh ows the location of the intergalactic 
absorption claimed by Nicastro et al. ( 20051 ) at z = 0.011. 
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Figure 8: The differential number of absorbers as a function of equi valent width 
forOVI A1035, OVII Ka, and OVIII Ka, bas ed on the model o f ICen fc Fang 
( 20061 ). Absorption measurements from OVI ( Danforth Shull . 20051 ) are in 
good agreement with the model. To detect the OVII and OVIII lines, one needs 
either path lengths of Az = 1 and a < 2 mA, or shorter path lengths (Az = 0.1) 
but a < 0.2 mA. This cannot be achieved with current instruments but should 
be achievable with future telescopes. 
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Figure 9: The soft excess emission can be expressed as the fractional amount of HI 
that one would need to add or remove to a spectral fit, assuming the ambient hot 
temperature for the cluster. Negative values indicate a soft excess. This quantity 
is strongly correlated with the Galactic 1/4 keV X-ray background, showing that 
the presence of an extra soft X-ray emiss ion component is due to in c orrec t removal 
of the Galactic soft X-ray background ( Bregman Llovd-Davies . 20061 ). 
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Figure 10: T ne dust extintion map of the Coma cluster region (from 
Schlegel et al.l fl998i ) . in units of magnitudes of extinction. A circle of 3 Mpc 
radius (1-7°) is cent ered on the Coma c luster . The crosses show the locations of 
the regions used by Finoguenov et al. ( 20031 ) to determine the properties of an 
X-ray filament. The Coma cluster happens to lie in a region of particularly low 
extinction, which may let through additional Galactic halo emission that could 
be mistaken for a feature intrinsic to Coma. 
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Figure 11: Two Chandra images of NGC 891 in the energy band 0.4-1.0 keV, 
where the region of high optical extinction is delineated by the red polygon. 
Each pixel is 0.5", so the field is 3.3' high and e ach dot is a single photo n. The 
X-ray shadow seen in the first observation (left, iBregman Irwin 2002 ) is not 
confirmed with more recent data (right). 



